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A new dynamics more natural than that proposed by Bouchaud and Dean is introduced to 
the Generahzed Random Energy Model, and the master equation for the dynamics is solved 
exactly to calculate the time correlation function C(i + t„, iw)- Although our results are very 
similar to those obtained by Bouchaud and Dean qualitatively, the exponents for power law 
relaxation are different. The Zero-Field-Cooled magnetization is also calculated with a relation 
between the correlation function and the response function which holds even if the relaxation is 
non-equilibrium. The validity of these analytic results are confirmed by numerical simulations. 
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§1. Introduction 

The aging, dynamical behavior largely depending on the history of system, is one of the most 
striking phenomena in the complex systems such as spin glasses, glasses, polymers and proteins, and 
many theoretical models have been proposed to explain these phenomena. The Generalized Random 
Energy Model is one of such models. Bouchaud and Dean calculated the correlation 

function C{t + t„,t„) exactly and showed that t/t^-scaling holds in this model^ although the 
dynamics they introduced seems to be a little unnatural. We recently introduce more natural 
dynamics to this model, and carried out simulationsi^ on aging phenomena with temperature 
variations. As a result, we revealed that the fundamental features reported in these complex 
systems, such as the memory effect and the rejuvenation, were well reproduced along the hierarchical 
picture discussed in ref. ^,[^,0. 

In this paper, we calculate the exact correlation function for dynamics used in ref. ^ and we 
confirm that while the t/tw-scaling holds as ref. ^ our result is different from theirs on exponents 
for power law relaxation. By introducing magnetizations to this model,@^ we also calculate the ZFC 
magnetization Mzfc(^; ^w)^^^ which is measured by quenching the system from above Tc to below 
and applying a weak magnetic field during t after a waiting time t„. We estimate Mzfc(^)^w) by 
using a relation between the response function R{t,t') and the correlation function C{t,t 

R{t,t') = -^^C{t,t'), (1.1) 
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where T is the temperature. This relation is vahd when assuming that the external field H is small 
enough and the transition probability from a state /? to a is modified by a factor exp(— ^^'^ ), where 
Mp is the magnetization of a state (3. This condition is satisfied in the present way of introducing 
magnetizations.^ Although this relation is very similar to the Fluctuation Dissipation Theorem 
(FDT) for equilibrium relaxation (actually it reduces to FDT if time homogeneity holds), it is valid 
even if the relaxation is non-equilibrium. These analytic results on the correlation function and the 
ZFC magnetization are checked by numerical simulations. 

The organization of this manuscript is following. In section 2, we introduce the GREM. In section 
3, we calculate the correlation function and the ZFC magnetization analytically and check these 
results by simulations. In section 4, we summarize this paper with some discussions. 

§2. Model 

The GREM is schematically shown in Fig. 1. The bottom points represent the accessible states 
of this system. This model is constituted by piling up L layers hierarchically. For simplicity, we 
hereafter consider the case that each branching point has branches. 

We say that d{a, (3), the distance between two states a and /3, is k if the lowest common ancestor 
of these two states is in the k-th layer counted from the bottom, e.g., in Fig. 1 d{a,f3) = 1 and 
d(a,7) = 2. It is assumed that overlaps between states depend only on the distance and they 
decrease with increasing distance. Therefore qk, the overlap between two states with d{a,(3) = k, 
satisfies the relation 

qo> qi> ■■■ > qL-l > qi = 0. (2.1) 

We hereafter denote as hk{Bk^i) the barrier height from a branching point in the k — 1-th layer, 
to its parent branching point. In the GREM, hk is given randomly and independently 
according to the distribution 

1 h 

Pk{hk)dhk = ——exp[-—^]dhk, (2.2) 
where Tc{k) is the transition temperature of the k-th layer and satisfies the inequality 

T,{1) < Te(2) < • • • < T,{L). (2.3) 

For dynamics, we employ Markoff process similar to that used in ref. ^. At first, the system 
is activated from a state /3 to Pk (the k-th. ancestor of /3) with probability W{f3, k) in unit time. 
After the activation, the system falls to one of all the states under /J^ with equal probability. We 
choose the uniform distribution for the initial condition to consider the situation that the system 
is quenched from an infinitely high temperature. 

The difference between the present dynamics and the one by Bouchaud and Dean exists in the 
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hopping probability W{P, k). Inspired by the model proposed by Yoshinojll' we give W{P, k) as 

1 



Wif3,k) ^ 



(2.4) 



where tq is a microscopic time scale, = exp(hk/T) and tl+i(/3l) = c«. Note that the first factor 
in the right hand represents the probability that the system can be activated from /3 to /3fc and the 
second one represents that the system can not be activated from to (3k+i- Hereafter tq is used 
as the unit of time and set to 1. From eq. ( |2.2[) , the distribution of can be written as 



Pk{fk)dfk = ^ , i dffc (ffc > 1), (2.5) 

where Xk = T/Tc{k). 

On the other hand, Bouchaud and Dean give the hopping probability as 

W{P,k)= ^ (2.6) 

T^On[Pk-l) 

But this hopping probability is not so realistic because, e.g., W{(3, A; + 1) > W{f3, k) if hk+i{f3k) < 
hkiPk^i)- 

Now let us introduce the following two functions. The first one ,Il{k, a, t), is the probability that 
the system which initially stays at a reaches a state with the distance k at t. This function is 
defined by the Green function G/^a{t), i.e., the probability that the system which initially stays at 
a reaches /? at t, as 

Il{k,a,t)= J2 Gf,a{t), (2.7) 

I3:d{a,l3)=k 

where the sum of this equation is taken over all the states (3 which satisfy the condition d{a, (3) = k. 
The next function, P{a, B,tvj), is the probability that we find the system at a state a at time 
when the initial condition is located anywhere under B, and is defined by the Green function as 

P(q, B, t^) = -i- ^ G^p{t^), (2.8) 

where N{B) is the number of states under B. Using these two functions, we can represent C{t + 
t„, tw), the correlation between t„ and t + t^, as 

) = E '?fcEn(^,«,i)^(«,^top,tw). (2.9) 

k=0 a 

We hereafter calculate the correlation function along the following strategy: 
(i) Calculate the Laplace transform of P[a,Btop,t^) 

POO 

P{a,Btop,E)= dt^exp[-Et^]P{a,Btap,t„). (2.10) 
Jo 
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(ii) Calculate the Green function and estimate 

/•OO 

fl(k,a,E')= dm{k,a,t)ex.p[-E't]. (2.11) 
Jo 

(iii) By taking the sum for a, calculate 

L-l 

C{E' + E,E) = Y,QkYl a, E')P{a, Btop,E). (2.12) 

k=0 a 

(iv) Carry out the inverse Laplace transformation of C{E' + E, E) to obtain C{t + t^, t^)- 

§3. Analytic Calculations and Simulations on the GREM 

3. 1 Derivation of the formal solution 

Now let us consider the event that the system initially stays at /? and reaches a at time tw The 
probability for this event is P{a,(3,tvi) = Gap{tw)- There are two possibilities for this event to 
happen: 

(i) a = (3 and the system is never activated during t^. 

(ii) The system is activated to (3^ at t' (< t^) and reaches a after that. 

In the case (ii), the probability that the system reaches a after the activation is P{a, I3k,t^ — t') 
because the next state is chosen randomly from all the states under Pk (recall the definition of 
P{a, B,t„)). Therefore, we obtain the following integral equation 

P{a,P,t^)=6^f,exp[-Jf-]+Y^ f'^ dt'W{(3,i)exp[-J—]P{a,P„t^-t'). (3.1) 

The Laplace transformation of this equation leads us to 

P(^. p ^a/3Ti(/3) I v ^(A^)n(/3) p, 



If we substitute the Laplace transformation of eq. ( p. 81) into eq. (p^), it becomes linear equa- 
tions of the Green functions, so that we can calculate the Green function P{a, /3, E) itself and 
P{a, Btop, E). The calculation is described in appendix and the results are 

P{a,B,,^,E) = rinn^ 7( V ^^'^^ 

and 

En{(3) + 1 ^^N^{En{a) + l}{Efm + l}ZkmUi=iZi(ai)Zi{f3iy 
where the function is defined as 

ZkiBk) = 1 - [1 - {fk+iiBk)}-^]Xk{Bk) {0<k< L), (3.5) 
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and Xk is defined by the following recursive equations, 

XoiBo) ^ (3.6a) 

3.2 Calculations of the correlation function 

Let us calculate Xj.{Bii) in the explicit form. From eq. (3.6), we can expect that Xj.{Bi^) does not 
depend on B^ in the limit — > oo (we hereafter simply denote as X^) and is calculated recursively 
as 

Xk = XkI{Vk,Xk) {l<k<L), (3.7) 

/(^,x)-r^^, (3.8) 

Ji u Vu + l 

where 

V, = (3.9) 

In the following calculations, we consider the case xl < • • • < a^i < 1- In the assumption E <^ 1 
(equivalent to iw » 1 ), we find 

Xk^l- CkE^'' il<k< L), (3.10) 



where 

k 

Ik = n^'' (3.11a) 

1=1 

{ r(i-xi)r(i + xi) {k = i), 

Ck = { (3.11b) 

[ r(i-xfc)r(i + xfe)C7,"^i {k>2). 

By using this result to eqs.(|3.5D and (|3.3|), P{a, Btop, E) is explicitly given as 



P(a, Stop, E) « ^ ^^n^r (a^_i) 

^ CL7V^{i?fi(«) + l}nf=inCii?^'fm(«0 + l] 
where we have used the fact f^+i = cx3. 

Next we estimate n(A:, a, ii/') defined by eqs. (|2.7| ) and ( |2 . 1 1]) . By taking the leading order of 

in eq. (^), we find 



{i?'n(a) + 1} nf=i Ziiai) . , {ii;'fi(/3) + 1} nti' ^/(A 



f3:d{a,l3)=k '^^t^' ' "Jll/=1 ^'VP^J 

fi(a)VF(a,A;) ''"^ 

«=1 



{S'ri(Q) + l}nr=i^KaO 



C.n(a)M.(a,.)i.-- (l<fe<L). (3.13b) 



{E'n{a) + l}Ui=iZiiai) 
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We are now in the position to calculate eq. ( |2.12 ). At first, let us define a function 

CkiE' + E,E)= J^Pia, Stop, E)U{k, a, E'). (3.14) 

a 

For A; = 0, this function is estimated as 

Co{E' + E,E)^ ^—xiJ{E, E'; xi - 2) J] xi+iI{CiE^' , xi+i - 1) 

1=1 

= 9{E',E;j,) + (3.15) 

where 

/■°° ds 

J{A,B;u)^ - (3.16) 
Ji s [As + l)[Bs + 1) 



and 



Similarly, for k > 1 



j^fa—lj^—a 

g{E',E;a)^ ^ _ ■ (3-17) 



17— 7i 1/~< k—1 

Ck{E' + E,E)^ ^-xiJ {E, E';xi-l)l[ xi+iJ {QE^^ , QE'^' ; xi+i - 1) 

xxfc+i I J (Cfc^;^'' , CfcS'^'= ; Xk+i - 2) - J (C7fc£;^^ , CkE'""^ ; rcfe+i - 1) | 

L-l 

X n X^+iI{CmE^^,X^+i-l) 
m=k+l 

^g{E',E;jk+i)-g{E',E;-fk), (3.18) 

where we have used the fact J {CkE""^ , CkE'"'f' ; x^+i - 1) is neghgible to J {CkE"^'^ , CkE'""' ; Xfc+i - 2). 

Now let us carry out the inverse Laplace transformation of eqs.( |3.15|) and ( |3.18| ). At first, ^^^]_^^ 
is transformed to 1. As for g{E' , E; a), the transformation from E' to t leads us to 

g{t, E; a) = ~^ f ds{t - s)-^E~^e''\ (3.19) 
r(l - a) Jo 

Then by transforming from S to t„, we find 

git, t^; a) = f (j^; «) " 1, (3-20) 

where 

fix; a) = —-^ d^(l - ur-'u"^. (3.21) 

1 (Q;)r(l - a) Jx 

From these results, we finally obtain 

L-l 



C(i + t^,tw) = X! IkCkii + tw, tw) 

k=0 



6 



The asymptotic behavior in regimes t <C tw and t ^ is 

1 - (t/tw)'~^^ {t « tw), 

Co(t + iw,tw) ~ <! (3.23a) 

(t/t^)-^' (t » tw), 

(t/tw)^"^^ (t«tw), 

Cfc(t + tw,tw) ~ { (3.23b) 

(iAw)~^"+^ (i»iw) (1<A:<L-1). 

As a result, 

go - (go - gi)(tAw)^"^' (t<tw), 
C(t + tw,tw) ~ <( (3.24) 

gL-i(tAw)"^^ (t»iw). 

3. 3 Calculations of the ZFC magnetization 

Now let us explain how to introduce magnetizations to the gremS We assign the value of 
magnetization to a state a as 

Ma = Moia) + Miiai) + ■■■ + ML-i{aL-i), (3.25) 

where Mkioik) is a contribution from the branching point a^, and is given from a distribution 
function T)k{M.k) as an independent random variable with zero mean. Note that A4n{n = k,k + 
1, • • •) are common between two states a and (3 if d{a, (3) = k. This means that there is a correlation 



of magnetizations between two states, whose averaged value MaMp is 



L-l 



M^Mp = J2 Ml {d{a,P) = k), (3.26) 



n=k 



where denotes the variance of T>n{Mn)- Therefore we can obtain the correlation function of 
magnetization by setting q^^ in eq. ( 3.22| ) to X]n=fc -^n- 



The effect of magnetic field H is considered by changing hi{f3) into /ii(/9) + HMp. As a result, 
the hopping rate W{(3, k) and the transition probability from (3 to any states are modified by a 
factor exp(— //Afg). Therefore, as mentioned in § we can use the relation eq. ( |1.1| ) to calculate 

MzFc{t,t^) as 

MzFc(i, t^) = H dt'R{t + t„, t') 

fe=o r(7fc)r(i - Ik) Jj^ 

3.4 Simulations 

To check the validity of analytic results obtained above, we carry out simulations on the GREM 
with L = 3, xi = 0.6, X2 = 0.4 and X3 = 0.3 ( see ref. |^ for the details of the simulation). We 
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first observe Ck{t + t^j^w)) the probability that the distance between the states at t„ and t + t„ 
is k, for = 10^, 10^, . . . , 10''. We take a random average over 10^ samples. In Fig. 2, we plot 
Ckit + twj^w) as a function of t/t^f with the analytic result shown in eq. ( |3.22| ). We can see that 
the data are rather consistent with the analytic result. 

Next we observe the ZFC magnetization MzFc(i; ^w) in the case that H/Tc{3) = 0.1 and 7W| = ^ 
{k = 0, 1, 2) for = 10^, 3 x 10^, 10^, . . . , 10^. The number of samples used for a random average 
is 10*^. In Fig. 3, the scaling plot of AfzFc(ii^w) is shown with the analytic result, eq. (3.27). The 
validity of the analytic result is well confirmed. 

§4. Summary and Discussions 

In this manuscript, we have introduced more realistic dynamics than that used in ref. |^ to the 
GREM, and calculated the time correlation function exactly. The results obtained by Bouchaud 
and Dean and ours are very similar in the sense that t/t^ scaling holds in both calculations and the 
asymptotic behavior shown in eqs. (3.23) and (|]2|) is the same. But in ref. |3|, all the exponents 
7fc (A; = 1, . . . , -L) in eqs. (3.23) and ( |3.24| ) are replaced with xj.. The reason for this difference on 
the exponents is that, as seen from eqs. (2^) and ( |2.6| ), the hopping rate W{/3, k) depends on fj(/?) 
{i = \, . . . ,k) in our calculation, but only on fk{l3) in ref. |^. 



We have also calculated the ZFC magnetization with the relation eq. (1.1) and shown the validity 
by simulation. From this result, we have confirmed that eq. ( [L.l| ) is really valid in the GREM. 
But, as mentioned in § |l], this relation holds only on the specific models and kinetics. Actually, 
Cugliandolo and Kurchan0 proposed a different relation between the correlation function and the 
response function on the SK model and confirmed it numerically. It is very interesting subject to 
investigate how these two functions are related in other systems in non-equilibrium. 

Now let us comment on two main reasons which seem to make this model solvable. The first one 
is that the number of branches in each layer is very large. For example, as seen from eqs. ( [j.3D , 
(^) and (3.6), P{a, Btop, E) depends on the full structure of the tree (f^ of all the branches, to 
be specific). But this quantity is self-averaging and only ffc(afe_i) {k = 1,2, . . . , L) dependence 
remains in the limit oo. The second one is the choice of the hopping rate W{l3,k). From 

eq. ( |3.12|) , we can see that P{a, Btop, E) depends on fA;(aA:_i) as the following form 

L 

P{a, Stop, ^) = n ^k{E, fk{ak-i)}. (4.1) 

k=l 

This means that P{a,Btop,E) depends on fi^{ak-i) independently. As seen above, this property 
makes our calculation simple. We confirmed that this property relies on the choice of W{(3, k) and 
does not hold if, for example, we set it as 

WiP,k)= ] (4.2) 

llj=i^Urt-i) 
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In this manuscript, we have succeeded to calculate the Green function of this model. By using it 
directly, we may be able to carry out simulations as done in ref. |^ in longer time regime. 
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Appendix: Derivation of eqs. (3.3) and (3.4) 

We first prove the following equations: 



P{a,Bk,E) 
b{Bk,i) 



[1 - {n+i[Ai{Bk)]}-^]XkiBk 



+ b{Bk,i)Pia,MBk),E) {0<k<L), (A-la) 

(A-lb) 



ZkiBk)UUk+iri[Ai^iiBk)] ' 

where ao = a. The notation Am{Bn) {m > n) stands for the branching point which is the ancestor 
of Bn and is in the m-th layer. The functions and are defined by eqs. (3.6) and (|3.5| ), and 
Yfc is defined as 

1 



E+l'' 



Yo{a) 



[l<k< L). 



(A-2a) 
(A-2b) 



We can easily check that eq. (A-1) is reduced to eq. ( p.2| ) for A; = 0. Next we can prove eq. (A-1 
in terms of mathematical induction by substituting eq. (AT) for k = n into the relation 

Pia,Bn+uE) = ^ P{a,Bn,E) (1 < n < L - 1). 



(A-3) 



B„£B„+i 



Equation. ( |3.3D is derived from eq. (AT) for k = L. 
For other branching points, P{a, Bf^, E) is given as 



P{a,Bk,E) =Y'^(.Bk,i)Yi{a)6A,(Bk),a, 

1 



i=k 



a{Bk,i) 



Zk{Bk) 

Xk{Bk)[l - {f^+i[Ai{Bk)]}-^] 



[ Zk (Bk) m=k+i Zi [Ai iBk)]fi (Bk)] 



{0<k< L), 

{i = k), 
{i > k). 



(A-4) 



(A-5) 



From these equations for = 0, we can derive eq. (|3.4|). We hereafter prove these equations 



recursively. For k = L, they are obviously valid from eq. ( |3.3|) . Next we show that eqs. ( A-4 ) 



and ( A-5 ) are valid for A; = n — 1 if they hold for k > n. By substituting eq. ( [A-4[) into the right 
hand of eq. (A-la) for k = n — 1, we find 

L 



where 



P{a,Bn-i,E) = ^" 1"" 1 +YS{Bn-i,i)Yi{a)6A,(B„^^),a, 



S{Bn^ui) = J2KBn~uk)a[AkiBn^i),i], 

k=n 



(A-6) 



(A-7) 
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Therefore eqs. (A^) and ( |A-5| ) are proved for /c = n — 1 if S{Bn-i,i) = a(i?„_i,i) (i > n). By 
substituting eqs. (A-lb) and ( A-5 ) into eq. ( A-7 ), it is shown that 



Zn-l{Bn-l) n m=n '^fn 



\Zi\MBn-l) 

a{Bn^i,i), 



i-l 

+ E 

k=n 



1 - Zk[AkiBn-l)] 

nu^/[^K^n-i)] 



(A- 



where we have used eq. 
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FIGURE CAPTIONS 



Fig. 1 The structure of the Generahzed Random Energy Model with L = 2 and N = 5. 

Fig. 2 Ck{t + tw)*w) versus t/t^ for t„ = 10^, 10'^, . . . , lO"^. The sohd hne indicates the analytic 
result obtained in eq. ( |3.22D . 

Fig. 3 MzFc(i,iw) versus t/t^ for t„ = 10^,3 x 10^, 10'', . . . , 10^ The solid line indicates the 
analytic result obtained in eq. ( |3.27 ). 
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